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We report on recent progress in the calculation of Bk using HYP-smeared improved staggered 
fermions on the MILC asqtad lattices. We have added measurements on fine (a ~ 0.09 fm) and 
superfine (a ~ 0.06 fm) ensembles at different values of the light sea quark mass (ami), as well as 
increased the statistics on some other ensembles. We find that the results on the fine lattices show 
a significantly stronger amg dependence than those on the superfine and coarse (a ~ 0.12 fm) 
lattices. We discuss different methods for accounting for these new results when doing the ami 
and continuum extrapolations. 
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1. Introduction 

The calculation of the kaon mixing matrix element Bk is one of the successes of lattice QCD. 
Results with all errors controlled are available with several fermion discretizations [JjJ The value 
for Bk is a key input into the standard model prediction for e K - At present there is an 3a tension 
between this prediction and the experimental value (if one uses the "exclusive" V c b obtained using 
lattice QCD) [|], ||, Thus it is important to continue to further improve the lattice calculations. 

In this proceedings we provide an update on our results for Bk- These are obtained using 
improved staggered fermions, specifically HYP-smeared valence quarks on asqtad sea-quarks. We 
describe here results obtained using chiral extrapolations based on SU(2) staggered chiral pertur- 
bation theory (SChPT), which is our most reliable procedure [[| Q]. Details of the fitting functions 
and the analysis method are given in Ref. [|L 0], and are not repeated here. 

Table [l] lists all the ensembles on which we have calculated Bk, and notes which results have 
changed in the last year. In particular, since Lattice 2011 we have accumulated much higher statis- 
tics on the F2 and SI ensembles (in addition to a small increase on the Ul ensemble) and added 
new measurements on ensembles F3, S2 and S3. Of these, the most important updates are those on 
ensembles F2, F3, S2, and S3, since they give information on the light sea quark mass dependence 
which was previously lacking on the fine and superfine ensembles. We focus on this dependence 
here. 



a (fm) 


ami /ani s 


geometry 


ID 


ens x meas 


status 


0.12 


0.03/0.05 


20 3 x 64 


CI 


564x9 


old 


0.12 


0.02/0.05 


20 3 x 64 


C2 


486x9 


old 


0.12 


0.01/0.05 


20 3 x 64 


C3 


671 x 9 


old 


0.12 


0.01/0.05 


28 3 x 64 


C3-2 


275 x 8 


old 


0.12 


0.007/0.05 


20 3 x 64 


C4 


651 x 10 


old 


0.12 


0.005/0.05 


24 3 x 64 


C5 


509x9 


old 


0.09 


0.0093/0.031 


28 3 x 96 


F3 


949x9 


new 


0.09 


0.0062/0.031 


28 3 x 96 


Fl 


995x9 


old 


0.09 


0.0031/0.031 


40 3 x 96 


F2 


959x9 


update 


0.06 


0.0072/0.018 


48 3 x 144 


S3 


593 x9 


new 


0.06 


0.0036/0.018 


48 3 x 144 


SI 


749x9 


update 


0.06 


0.0025/0.018 


56 3 x 144 


S2 


799x9 


new 


0.045 


0.0028/0.014 


64 3 x 192 


Ul 


747 x 1 


update 



Table 1: MILC asqtad ensembles used to calculate Bk- ami an d am s are the masses, in lattice units, of the 
light and strange sea quarks, respectively, "ens" indicates the number of configurations on which "meas" 
measurements are made. Note that the numbering of the ID tags on the fine and superfine lattices do not 
follow the ordering of ami- 



2. Chiral fits 

In our numerical study, our lattice kaons are composed of valence (anti)quarks with masses 
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Figure 1: X-fits to B^(l /a) for the F3 and F2 ensembles. The red diamond is explained in the text. 



m x and m y . These are, respectively, the masses of the valence d and s quarks. On each MILC 
ensemble, we use 10 valence masses: 

am x ,am y = am s xn/ 10 with n = 1,2,3, 10, (2.1) 

where am s is the nominal strange sea quark mass. In our standard fits we extrapolate to am^ hys using 
the lowest 4 values for am x (the "X-fit" — done at fixed am y ), and then extrapolate to /Hj hys using 
the highest 3 values of am y ("Y-fit"). As described in Ref. [0], these choices keep us in the regime 
where we expect next-to-leading order (NLO) SU(2) ChPT to be reasonably accurate. The X-fits 
described here are done to the form predicted by NLO partially quenched SChPT (which is given 
in Refs. [||, |8|]), augmented by analytic and generic non-analytic terms of next-to-next-to-leading 
order (NNLO) and a single analytic term of next-to-next-to-next-to-leading order (NNNLO). We 
use the Bayesian method to constrain the higher order coefficients, as described in Refs. f|, 
The am y dependence (which is not controlled by ChPT) is very close to linear and we use a linear 
fit for our central values. We dub this entire fitting procedure the "4X3Y-NNNLO fit", and use it 
for our central values. 

The X-fits on the new fine ensembles are shown in Fig. |l|, while those on the new superfine 
ensembles are shown in Fig. ^. In each panel, a fit to the SU(2) SChPT form is shown, with Xp the 
squared-mass of the Goldstone-taste valence xx pion, whose mass is very close to linear to am x . 
Also shown (as the red diamond) is the result obtained after (a) extrapolating am x — > fl/n' hys , (b) 
removing the known taste-breaking in the pion masses appearing in the chiral logarithms, and (c) 
setting the light sea-quark mass to its physical value in the chiral logarithmic terms. See Ref. [||] 
for details of this procedure. We have incorporated finite volume corrections, as predicted by NLO 



SChPT, into the fitting function [|T0p. 

A notable feature of both figures is that the difference between the extrapolated/corrected result 
and the fit curve is larger on the ensembles with larger ani( (F3 and S3). This is due almost entirely 



3 



Recent progress in staggered Bk 



Weonjong Lee 



0.52 



0.5 



m 



0.48 



0.46 




4X3Y-NNNL0, X-fit 



_L 



0.52 



0.5 



0.48 



0.46 




4X3Y-NNNL0, X-fit 



0.05 0.1 0.15 0.2 



X p (GeV^) 
(a) S3 



X p (Gey**) 
(b) S2 



Figure 2: X-fits to Bk{1 /a) for the S3 and S2 ensembles. Notation as in Fig. |T| 
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Figure 3: Bk(p = 2GeV) vs. Lp (GeV 2 ), where Lp is a squared mass of the Goldstone-taste pion composed 
of two light sea quarks (It). 



to the larger shift required to bring the sea-quark pion mass to its physical value on these lattices 
(correction (c) above). The other aspects of the extrapolation/correction are similar for the pairs of 
lattices. We return to this point below. 



3. Dependence On Light Sea Quark Mass 

In Fig. |3| we show the dependence of Bk on the light sea quark mass for all lattice spacings, 
including the Ul point for completeness. We use Lp as a proxy for ami, since these quantities are 
proportional to very good accuracy. We have run Bk to a common renormalization scale /I = 2 GeV 



4 



Recent progress in staggered Bk 



Weonjong Lee 



(in the MS scheme), so that the results are directly comparable. We expect, at NLO, only a linear 
behavior on hp. 

As can be seen from the plot, a linear fit works well for each of the three lattice spacings for 
which we have multiple values of ami. The results on the coarse lattices, which were presented 
in Ref. [Q], show a mild dependence on hp. The slope, —1/(2.9 GeV) 2 , is consistent with naive 
dimensional analysis, which predicts a magnitude of 0{\ /A 2 ) with A x Ft 1 GeV. The slope on the 
superfine lattices, which is a new result, has a slightly larger magnitude than that on the coarse lat- 
tices but is comparable. What is striking, however, is the slope on the fine lattices, whose magnitude 
is about twice that on the coarse lattices (slope w — 1/(2 GeV) 2 ). Its value per se is reasonable — 
what is unexpected is that it should differ so much from that at the other lattice spacings. Since we 
are in the regime where corrections are small (< 10%) in both chiral and continuum extrapolations, 
one would expect that the slope would depend only weakly on a 2 , and any dependence should be 
linear. Instead, our results appear to require a contribution to Bk proportional to ani((a 2 ) 2 , since 
only with a quadratic dependence on a 2 can the magnitude of the slope increase and then decrease 
as a 2 is raised. 

This peculiar behavior has only recently come to light and we are in the process of investigating 
it. One avenue we are following is to do a global continuum-chiral extrapolation to all results. This 
would incorporate the constraint that fit parameters in the chiral fits should depend smoothly on 
a 2 . Another approach is to remove the step (c) described above from our correction procedure and 
leave this to be dealt with by the ami extrapolation. This is possible because, at NLO in SChPT, hp 
does not appear in the argument of the chiral logarithms, but only in the prefactor. Pending a more 
complete understanding, we have considered the continuum extrapolation of results obtained after 
the extrapolation to physical ani(. 

4. Continuum Extrapolation 

In Fig. Q we show the results obtained on the coarse, fine and superfine lattices after the 
linear extrapolation to dmj hys plotted vs. a 2 . We also include the Ul result which, however, is not 
extrapolated to physical ami, since there is only one value of ami available for the ultrafine lattices. 
Thus it cannot be consistently used in the continuum extrapolation. 

Previously, when the high-statistics data from ensembles F2, F3, S2 and S3 were not available, 
we made the continuum extrapolation using ensembles C3, Fl, SI and Ul, all of which have 
close to the same sea-quark masses. We then included systematic errors accounting for the use of 
unphysical sea-quark masses (the errors being 1.5% and 1.3% for the light and strange sea-quark 
masses, respectively [|J]). We found a continuum behavior which was far from linear in a 2 (as can 
be seen from Fig. || — the fit corresponds to hp « 0.1 GeV ), and we could not obtain reasonable 
fits including the coarse ensemble. Dropping this, and trying various fits (along the lines described 
below) we estimated an extrapolation error of 1.9%. 

Using the data extrapolated to physical ami leads to a much smoother continuum extrapolation, 
as can be seen in Fig. |[ We can now obtain good fits including the coarse lattice result. As noted 
above, we can consistently fit only to the larger three values of a 2 . Following Ref. [Qj, we fit both 
to a linear dependence and to a five-parameter form containing dependence that we know enters 
due to incomplete operator matching and higher-order contributions. The explicit forms are given 
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Figure 4: Bk(}1 — 2GeV) vs. a 2 (100 x fm 2 ). See text for description of data points. Note that the Ul point 
(smallest a 2 ) has not been extrapolated to physical ami. Red, brown and blue curves and points show the 
lin3, lin4 and g4-a2g2-a4 fits, respectively. 

in Table [2], For the five-parameter fit we use Bayesian constraints for C2 — c$, with central values 
set to zero and spreads taken as 02 = (73 = CJ5 = 2 A 2 and 04 = 2, with A = 300 MeV. These fits are 
shown in Fig. |j; values for # 2 /d.o.f. and extrapolated Bk are listed in the Table. Also shown, for 
completeness, is a linear fit including the Ul point, labeled lin4. 



fit type 


#data 


fit function 


£ 2 /d.o.f. 


5*(2GeV) 


lin3 


3 


c\ +c 2 a 2 


0.0307 


0.5346(51) 


g4-a2g2-a4 


3 


c\ +C2C1 2 + C3a 2 a s (a) 










+c/s r a 2 (a) +csa A 


0.0298 


0.5178(133) 


lin4 


4 


c\ + cicr 


0.703 


0.5314(43) 



Table 2: Details of continuum fits. For the g4-a2g2-a4 fit with Bayesian constraints, the augmented 
^ 2 /d.o.f. is shown. Only statistical errors are shown for Bk- 

At this stage we are not ready to quote an updated result for Bk, since we clearly need to 
improve our understanding of the ami dependence. Nevertheless, we note that, taking the linear 
continuum extrapolation (fit lin3) the final value for Bk is 0.732 with a ~ 1% statistical error. 1 
This is, in fact, consistent with our previous result of 0.727 ± 0.038 [Q] (where the error is domi- 
nated by systematics). 2 Setting aside the matching error of 4.4%, we note that the combined am/ 
and continuum extrapolation systematics were estimated at 2.4%. This is certainly large enough 
to accommodate any small shift we might find from our new analysis. Thus we do not expect a 
significant change in our final result. We do hope, however, that the ami extrapolation error will 

'We have also taken the opportunity to update the value of r\ to 0.3117fm, as given in Ref. flll|]. 

2 The smallness of the shift in the central value can be understood from Fig. ^. Extrapolating at physical Lp there 
is a larger a 2 dependence than at Lp ~ 0. 1 GeV 2 (ignoring the coarse lattice point in the latter case). This difference 
approximately cancels the increase in the values one obtains when extrapolating to physical Lp. 
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come down. The total error will, at present, remain dominated by the matching error (due to our 
use of 1-loop perturbative matching). 

In summary, results at different light sea-quark masses on the fine and superfine MILC asqtad 
ensembles have uncovered an ami dependence that is hard to understand. An optimistic interpre- 
tation of the new results is that the difficulties we had previously when extrapolating at unphysical 
ami were due to artifacts introduced by our X-fitting and extrapolation procedure, while the new 
continuum fits (shown in Fig. |]) do not contain these artifacts. Analysis is underway to test this 
interpretation, and final results should be available soon. 
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